039893 


.Bi^williiiu^  iitit'  iuiiti^gii^ 


FINAL  REPORT 


Aagittt  1975  to  Jvify  1976 


peeaAiiM>r  1976 


AiMnriQKadr^^  Disirmtiaii  Unlimitod 


IhreptJridijt^  5 ■ 

'V-'?'.:'.  ■■  ■?■.  ■'.  '■.C'K^  ■■  ■ ■■■■'V’l^' ■^■■■■^  ■'■■  ■■  Cv..' . ■ 


^.11.  Aiil  tOROE  dCOE^ra 

^ v;5t.-  .•  r'--  V.  -cl''-'' '.  . 

Am^FOftdE^SYSWW  ; , . 


Slnii^xo  stM^  AIR 


EBRfORD*  i^UERApRdSETTB  6t7» 


'9ST' 


’•■  > - ■■•  - ■ * ‘ " :•■■•-.■  ./'■  ■ " •■ 


.S- 


;.-h 


■■?, 


f'.  '.  • 


...  t- 


rp-?uM-. 


!^;''vf''*f 


m 


'M 


M 


'^■vs  .-"^v-', 


■■y- 


ypm 


■i'^'S' 


WjlTS- 





.i»y 


yi,wr« 


1^ 

V|5  ,’'i, 

SaU".:i<;"'-'' . 

'■^'.fr.V^': 


';?T^  /it-. 


* ,5%  ^’’'J.. 


'J^'-Vv--.' 


‘^*SI 


.'/Jf-  ;.--..’St.V 
■ ^'' ■■■  i'-.J'  ..  ■■ 

. i - 


• ' J ; Th  .-J-  r" 

f •-.  ■.^•''■v  • 

. ,;-.>  • I 4^ 


f'.fp'' 


i'  ;[ 


■SiX 


, ^tQualtfied  rei|ue«toni  may  oteUin  addltlo^ 

OocumastaSm  Canter.  AU  others  ahould  apfdy  to  the  National  ^ 
Tec^ioieal  information  Service. 


St->' 


security  classification  of  this  page  rwhwi  Data  EntaradD 


PORT  DOCUMENTATION  PAGE 


READ  INSTRUCTI(»(S 
BEFORE  COMPLETING  FORM 


ENT'S  CATALOG  NUMBER 


JTUDY  TO  JffiVBLOT  GRADIOMETER  jJpMPENSATK»r  | 
^HNIQOE?  0 ~ ! 


Final 

Angmt  ^75— Ju 


Stanford  University  . 

Guidance  fc  Control  Laboratory 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Geophysics  Laboratory 
Hanscom  AFB,  Hass  01731 


4.  monitoring  agency  name  a AOORlSVif  dlllamtl  froai  Controlling  OllleoJ  I IS.  SECURITY  CLA 


Unclassified 


ISa.  DECLASSIFICATION/OOWNCRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  fol  IMo  ftapori; 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  fof  Ula  abatraet  aniarad  In  Block  20,  II  Mflotont  horn  RaporO 


aisfliflimv 


ft.  KEY  WOKDS  (C«9%tinu9  on  rorormo  old*  U noooooorf  mnd  Idontify  ty  Mock  mmbor) 


gravity  g;radiometer ; Hughes ^£fav?fy  Gradlometer  (RGG);  Bell  Aerospace, 

C. Stark  Draper  Lab.,  Evtvos  Units,  spin  axes,  sensors;  stress  relaxation, 
creep. 


10.  ABSTRACT  (Conthnio  on  raaaraa  alda  II  noeoooorr  and  Idontllr  Of  Mock  nuaibarj 


see  other  side 


DO  ,:2:**7,  1473  EDITION  OF  I NOV  OS  IS  OBSOLETE 


S/N  0102  LF  0144601 


ra 


security  CLASSIFICATION  OF  THIS  PAGE  fSfcan  I 


SeCUNITV  CCASSIPICATION  Or  THIS  HAGEOntMi  Om«  Enitnd) 


There  ere  three  eovlng  base  gravity  gradiometers  currently  under  devel- 
opaent.  The  Instruments  are  being  develed  at  the  Hughes  Research  Labs  ^Ref . TJ", 
the  Bell  Aerospace  Corp.  and  the  Charles  Stark  Draper  Lab,  The 

design  goal  for  each  of  the  sensors  is  1 EUtvBs  Unit  (EU).  The  group  includes 
sensors  designed  specifically  to  measure  the  gravity  gradient,  as  well  as 
sensors  which  utilize  existing  accelerometers  to  provide  a gradient  estimate. 
The  Hughes  and  Bell  instruments  rotate,  modulating  the  information.  This  ro- 
tation transfers  the  gravity  gradient  signal  to  a higher  frequency,  quieter, 
portion  of  the  spectrum,  and  can  separate  the  signal  from  some  sources  of  in^ 
strument  bias.  The  Draper  Lab.  sensor  measures  the  gradient  signal  at  zero 
frequency  and  uses  a sophisticated  flotation  suspension  system  to  isolate  the 
sensing  element  from  errors  induced,  by  rotation  and  Jitter. 

A system  of  at  least  three  instruments  of  any  one  type  is  required  in 
order  to  provide  a complete  gravity  gradient  tensor  estimate.  The  primary 
objective  of  this  paper  is  to  define  the  optimal  relative  orientations  of  the 
three  sensors  which  comprise  a minimum  gravity  gradient  measurement  system. 

In  addition,  the  effects  of  different  relative  sensor  orientations  when  a 
fourth,  redundant  sensor  is  included,  are  considered. 
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A portion  of  the  work  done  under  this  contract  has  been  combined 
into  the  Final  Report  of  the  previous  contract  In  order  to  keep  the 
material  on  model  fitting  as  an  Integrated  report.  That  report  repre- 
sents the  principal  parts  of  the  doctoral  thesis  of  E.  J.  Pelka.  Some 
work  that  has  been  done  since  that  thesis  study  Is  incorporated  as 
Chapter  I of  this  Final  Report,  which  represents  the  last  part  of  the 
studies  done  In  this  subject  area  under  the  current  contract  though 
not  all  of  it. 

The  literature  search  on  stress  relaxation/creep  reported  in  Chap- 
ter II  was  largely  negative.  Discussions  with  people  in  the  Materials 
Science  Department  at  Stanford  and  elsewhere  indicate  little  Is  known 
about  the  mechanism  and  nature  of  mlcrostralns.  Professor  John  Shyne 
has  recently  spent  several  months  In  Japan  interacting  with  a group 
doing  research  In  this  area.  It  Is  hoped  that  a better  understanding 
of  this  phenomena  will  Improve  our  ability  to  Interpret  long  term  var- 
iations In  gradlometers  which  have  been  observed  in  some  of  the  test 
data  reported  by  the  developers. 

Research  that  was  started  In  the  evaluation  of  small  accelerations 
to  be  used  for  the  on-line  correction  of  kinematically  Induced  errors  began 
with  disappointing  performance ; thus  it  was  felt  it  would  be  better 
omitted  In  this  Final  Report, 
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lOTRODUCTlON 

There  are  three  moving  base  gravity  gradlometers  currently  under  develop- 
ment. The  Instruments  are  being  developed  at  Hughes  Research  Labs  [Ref.  1], 
the  Bell  Aerospace  Corp.  [2],  and  the  Charles  Stark  Draper  Lab.  [3,  4].  The 
design  goal  for  each  of  the  sensors  is  1 Ebtvbs  Unit  (EU).  The  group  includes 
sensors  designed  specifically  to  measure  the  gravity  gradient,  as  well  as 
sensors  which  utilize  existing  accelerometers  to  provide  a gradient  estimate. 
The  Hughes  and  Bell  instruments  rotate,  modulating  the  information.  This  ro- 
tation transfers  the  gravity  gradient  signal  to  a higher  frequency,  quieter, 
portion  of  the  spectrum,  and  can  separate  the  signal  from  some  sources  of  in- 
strument bias.  The  Draper  Lab.  sensor  measures  the  gradient  signal  at  zero 
frequency  and  uses  a sophisticated  flotation  suspension  system  to  isolate  the 
sensing  element  from  errors  Induced  by  rotation  and  jitter. 

A system  of  at  least  three  instruments  of  any  one  type  is  required  in 
order  to  provide  a complete  gravity  gradient  tensor  estimate.  The  primary 
objective  of  this  paper  is  to  define  the  optimal  relative  orientations  of  the 
three  sensors  which  comprise  a minimum  gravity  gradient  measurement  system. 

In  addition,  the  effects  of  different  relative  sensor  orientations  when  a 
fourth,  redundant  sensor  is  Included,  are  considered. 


THE  GRAVITATIONAL  FIELD  AND  ITS  GRADIENT 

A gravitational  field  is  conservative  and  may  be  expressed  as  a 
scalar  function  cp  of  position  relative  to  the  attracting  body. 

The  gravitational  acceleration  "g  due  to  the  body,  at  any  point  "r,  may 
be  written  as 


g » '7(p(r)  , 


(1) 


! ^ 


Expressing  g in  a set  of  Cartesian  coordinates,  (1)  may  be  rewritten  as 


-1- 


k 


r-»-lT  fScD  S® 

" [sx  ay  azj 


The  spatial  gradient  of  the  gravity  vector  "g,  written  in  Cartesian 
coordinates  is  the  second  order  tensor 
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In  order  to  determine  completely  the  gravity  gradient  tensor  of  (3), 
all  of  the  independent  components  must  be  measured.  For  a conservative 
gravitation  force:  [5] 

^ X g = 0 ; (4) 

end 

V . g a 0 . (5) 


Equation  (4)  demonstrates  the  symmetry  of  the  gradient  tensor.  Equa- 
tion (5)  is  Laplace's  equation  and  requires  that  the  sum  of  the  diag- 
onal components  of  the  gradient  tensor  is  zero.  Since  the  gradient  tensor 
is  symmetric,  no  more  than  si.x  of  its  nine  components  are  independent. 

The  scalar  Laplace  equation  (5)  further  reduces  the  total  number 
of  independent  elements  by  one.  Hence,  only  five  of  the  nine  elements 
of  the  gravity  gradient  tensor  are  independent.  Complete  determination 
of  the  gravity  gradient  tensor  at  any  point  requires  only  five  measure- 
ments. 

To  follow  convention,  the  gradient  tensor  will  be  represented  by  f, 
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B. 1 Simplified  Gradlometer  Signal  Development.  Steady  State  Operation 


Approximation. 

Consider  the  ideal  dumbbell  gradlometer  oriented  with  its  spin  axis 

^ ~ y 

z horizontal  as  shown  in  Fig.  2.  The  vector  g is  the  gravity  vector 
s 0 

at  the  origin  of  the  sensor  axis  reference  system  and  defines  the  local 
horizontal  plane. 

At  any  of  the  masses,  m^^,  for  example,  the  gravity  vector  can  be 
written  in  terms  of  g^^  as  in  (7); 


*1,  ‘ 'o,  * aT*- 


J = x,y,z 


k = x,y,z 


Using  (6)  and  the  definition 

At  i (Ax,  Ay,  , (t 

the  gravity  at  m^^  may  be  written  as 

^ ^^1  • C9 

Since  At^  refers  to  the  displacement  of  m^^  from  the  origin  of  the 
sensor  reference  system. 


= Tj  , 


and  (9)  may  be  v/ritten  as 


-»  -*  =4  _* 

'i  ■ «o  * • 'l 


The  torque  % Induc’ed,  by  the  eravltatlonal  field  influence  upon  the 
odd  arm  may  be  written  as 


where 


Since  for  an  ideal  instrument  all  four  masses  are  equal,  and  since 


^3  = -^1  ' 


(12)  can  be  rewritten  as 


To  = 2mrj^  X {r  • 


(T  - Y ) • 7.  = 2mr^|(r  -r  )sin  29  + 2F  cos  29>.  (19) 

' o e'  s I yy  xx'  xy  J 


7,  = 


An  equation  similar  to  (16)  can  be  written  for  even  arm  induced  torques 
and  therefore  the  differential  output  torque  induced  by  the  gravity 
gradient  is 


T - T = 2m 
o e 


Along  the  sensitive  axis,  z^,  the  component  of  the  differential  torque 
vector  of  (17)  is 

^ (T  - T ) • z = 2mr^(2(r  - V )sin0  cos9 

' o e'  s \ yy  XX 

+ 2r^y(cos^9  - sln^0)|  (18) 
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By  separating  the  Inphase  (sin  26)  and  quadrature  (cos  26)  components  of 
output  Eq.  (19),  measurements  containing  three  of  the  five  Independent 
terms  of  the  gravity  gradient  tensor  are  available. 


£.  BELL  AEROSPACE  CORP.  GRADIO^!^ETER 

The  sensitive  elements  of  the  Bell  Aerospace  Corp. 's  gravity  gradio- 
meter  are  four  modified  Bell  Model  VII  pendulous  accelerometers.  The 
accelerometers  are  mounted  symmetrically  on  a spin  table  with  their  sensi- 

/N 

tive  axes  in  the  tangential  direction  defined  by  s,  where 


s 


w X 


(20) 


The  nominal  spin  table  rate  is  0.5  Hz.  A conceptual  diagram  of  the  Bell 
instrument  was  shown  in  Fig.  3,  The  instrument  senses  gravity 
gradient  in  the  X-Y  plane  according  to  the  equation 


) 


j 


[(3i  + a^)  - (ag  + a^)]  = 2d{(ryy-r^^)  sin(2wt) 

+ 2?  cos(2u)t)}  . (21) 

xy 


Modulation  of  the  output  signal  yields  an  ihphase  measurement  of  the  off- 

diagonal  gravity  gradient  tensor  component,  P , and  a quadrature  meas- 

xy 

urement  of  the  difference  between  two  of  the  diagonal  terms.  P - P . 

XX  yy 

Comparison  of  (19)  and  (21)  shows  that  the  output  equations  for  both 
of  the  rotating  gravity  gradiometers  are  equivalent. 


£.  DRAPER  LABS  GRAVITY  GRADIOMETER 

The  Draper  Lab's  instrument  senses  the  gravity  gradient  signal  at 
zero  frequency,  i.e.,  the  instrument  does  not  rotate.  A conceptual  draw- 
ing of  the  Draper  Lab's  spherical  gradiometer  is  shown  in  Fig.  4. 

The  spherical  gradiometer  uses  a spherical  float  having  dense  masses 
attached  at  two  points  which  provide  a dumbbell  mass  distribution  as  shown 
in  the  figure.  External  torques  are  sensed  about  each  of  the  two  perpen- 
dicular axes  (i,  J)  which  define  a plane  normal  to  the  between-the-weights 
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FIG.  4 SIMPLIFIED  REPRESENTATION  OF  DRAPER  LAB'S 

SPHERICAL  GRAVITY  GRADI(»IETER ; OUTPUT  AXES 
i AND  J. 


axis  (k).  The  flotation  fluid  is  inviscid,  providing  a degree  of  isolation 
from  small  amplitude  (Jitter)  rotations  of  the  case.  The  output  vector 
equation  for  the  spherical  g .''ometer  is:  [6] 


MEASURING  THE  GRAVITY  GRADIENT  WITH  ROTATING 
GRAVITY  GRADIOMETERS'*' 


In  Fig.  5,  it  is  assumed  chat  the  sensors  are  not  spinning  and  that 

and  that  the  sensor  axes  (x  , y , z ) are  respectively  coincident  with 

«i  J 

platform  axes  (x^,  y , z ).  The  separate  sensor  reference  system  orienta- 
tions can  then  be  arbitrarily  varied  using  a {1,2,  3}  Euler  rotation 


Since  the  signal  equations  are  the  same  for  both  the  Hughes  and  Bell 
sensors,  it  is  understood  that  the  analysis  of  signal  estimation  for 
the  Hughes'  sensor  applies  directly  to  the  Bell  sensor  also. 


FIG.  5 HUGHES'  SENSOR  SYSTEM  INITIAL  (MlIENTATION 

FOR  PURPOSES  OF  SPIN  AXIS  ALIGNMENT  SPECIFICATION 


sequence  defined  by  angles  {Oj,  Pj,  Yj}»  Subsequent  to  these  alignment 
rotations,  the  sensors  are  spun  up  and  measurements  taken.  The  output  of  each 
sensor  corresponds  to  (19),  with  each. censor  measuring  components  of  the 
gravity  gradient  tensor  expressed  relative  to  the  newly  specified  body 
(sensor)  reference  directions.  In  order  for  these  measurements  to  provide 
a complete  estimate  of  the  gravity  gradient  tensor,  all  six  outpqts  from 
the  three  sensors  must  be  related  to  a single  reference  frame.  For  con- 
venience, this  frame  has  been  defined  to  be  the  platform  reference  frame. 

It  is  easy  to  show  that 


^b/p^p’^b/p 


(23) 


Sensor  inphase  and  quadrature  signals  are  averaged  over  a 10  sec  sample 

period  to  provide  measurements  of  fP  - T and  [P  1 . as  speci- 

f yy  XX  ■'BJ  '•  xy*Bj 

fied  by  (24)  and  (25), 


i 


= 4mr  [i(P  -P  )]„, 

L IpJj  '-2'  yy  XX  ^BJ 


[A,Ij  = . 


(24) 

(25) 


y 


Using  (23),  the  right  sides  of  (24)  and  (25)  can  be  expressed  in  terms  of 
the  platform  referenced  gradient  tensor  for  arbitrary  Euler  angle  rota- 
tions {a,  P,  y}.*  Equations  (24)  and  (25)  can  be  rewritten  for  all  three 
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[W]  [r] 


The  6x6  matrix  M'  of  (26)  and  (27)  is  singular  due  to  the  Laplace 

relationship  of  r , p , r as  shown  in  (5),  The  effect  of  this 

xx'  yy  zz 

singularity  can  be  removed  either  by  eliminating  one  of  the  gradient 
tensor  principal  diagonal  terms  using  (5)  or  by  adding  the  Laplace  equa- 
tion to  (26)  as  a 7th  measurement  having  zero  error.  The  addition  of 
Laplace's  equation  as  a measurement  allows  all  of  the  terms  of  the 
gradient  tensor  to  be  treated  in  the  same  way  mathematically.  Combina- 
tion of  (5)  with  (27)  yields 


6X6 

1 1 *1  0 0 0 


[A]  = [P]  . 


[B]  i [m"^] 


the  least  squares  error  solution  estimate  of  [P]  is 


[f]  = BA, 


where 


A^  ^ A . V . 


Matrix  V is  a 7X1  vector  representing  measurement  noise.  The  first 

six  terms  of  vector  V are  taken  to  be  white,  zero  mean  Gaussian  random 

2 

variables  having  variance  1 (EU)  . The  7th  component,  V^,  is  exactly 
zero.  Furthermore, 


«(V  V‘)  « 


°1X6  ' 


7X1  • 
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Equation  (33)  states  that  all  six  measurements  have  uncorrelated  errors 
of  variance  1 EU^.  With 


where 


A ^ 

r = r “ r . 


^X6  * ‘t!*  *'1  ■ 


P = B R B 


R = S[V  V ] . 


Given  (36),  the  problem  is  to  find  angles  (a,  y)j  (which  affect 
B)  so  that  a scalar  cost  function  J,  which  is  related  in  some  way  to 
covariance  matrix  P is  minimized.  One  reasonable  specification  for  the 
cost  function  J is  [cost  function  suggested  by  Prof.  John  V.  Breakwell] 


K 


J = Tr{s[r  r ]} 


where  F is  defined  by  the  3x3  matrix  composed  of  the  errors  in  the 
estimate  of  With  f defined  by  (34), 


''l  ''4  ''S 


''2  ''e 


^ ^ ''a 


M ^ ^ ^ fir 
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The  off-diagonal  terms  of  (40)  are  all  zero,  and 


Tr  = Jj  rj  + 2 Pj  . 


T = orthonormal  transformation  matrix  relating  the  expression  of  a 
vector  in  reference  system  A to  its  expression  in  reference 
system  B, 


and  since 


Using  (42)  it  is  easy  to  show  that 


» fiS  T 

Pa  = TTgT 


Hence , 


= if 

A A B 


and  since  the  trace  is  invariant  under  similarity  transformation, 
TrCS[r.r,]}  = Tr{s[r  r_]}  . 


Therefore  the  cost  function  J defined  by  (38)  is  invariant  under  linear 
transformation.  This  invariance  allows  the  cost  function  tp  provide  a 
scalar  measure  of  the  quality  of  the  gradient  estimate  without  regard 
for  the  coordinate  system  to  which  the  measurements  are  referenced. 

The  system  shown  in  Fig.  6 is  taken  as  the*  baseline  gradient  measure- 
ment case.  The  spin  axis  of  each  sensor  is  parallel  to  a different  plat- 
form reference  direction.  For  this  sensor  system  orientation,  J had 

2 

the  value  8.667  EU  , The  Euler  angles  corresponding  to  this  set  of 
sensor  orientations  are  listed  in  (46), 
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FIG.  6 BASELINE  HUGHES'  GRADICMETER  GRADIENT  MEASURING  SYSTEM 
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This  particular  set  of  sensor  spin  axis  orientations  minimized  J for 
all  sets  of  sensor  spin  axis  directions  considered. 

A test  of  Hughes'  sensor  system  gravity  gradient  estimation  for  the 
case  where  the  spin  axes  were  not  orthogonal  is  depicted  in  Fig.  7, 


FIG,  7 CLOSING  THE  SPIN  AXIS  UMBRELLA  ON  x 


This  test  was  initialized  by  symmetrically  placing  the  three  sensor  spin 

A»  /S  /S 

axes  about  x in  the  y — z plane.  The  sensor  spin  axes  were  then 
P P R ^ /s 

uniformly  rotated  out  of  the  y - z plane  towards  x through  an  in- 

P P P 

creasing  angle  ^ in  a manner  similar  to  closing  an  umbrella.  The 
results  of  this  test  are  given  in  Fig.  8.  As  the  umbrella  angle  increased 


from  5 deg  to  35  deg,  the  cost  function  decreased  from  13 


.9  EU^ 


to  its 


minimum  value  of  8.667  EU  at  35.26  deg.  At  that  angle,  the  spin 
axes  were  orthogonal.  As  the  umbrella  continued  to  close^  the  cost  func- 
tion again  began  to  increase,  principally  due  to  a loss  in  the  accuracy 
of  the  r estimate. 

XX 
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The  cost  function  analysis  was  extended  to  ibe  4-lnstrument  case.  A 
fourth  instrument  provides  protection  against  loss  of  mission  due  to  the 
failure  of  a single  gradiometer.  Several  four-instrument  operating  con- 
figurations were  considered.  These  were: 

(1)  All  four  spin  axes  w’ere  symmetrically  placed  about  x in  the 

A.  ^ 

y - z plane.  The  umbrella  was  closed  with  all  four  spin 
P P 

axes  collapsing  towards  x^; 

(2)  One  spin  axis  was  oriented  along  x^.  The  three  other  spin 

axes  were  placed  symmetrically  about  x in  the  y - z 

P P P 

plane.  The  umbrella  was  closed  as  in  (1}. 

The  results  of  these  tests  appear  in  Figures  9 and  10.  Figure  9 shows 

that  for  the  sensor  configuration  defined  by  case  (1),  an  umbrella  angle 

of  35  deg  is  optimal  for  both  the  full  four-sensor  case,  the  cost  func- 

of  one  instrument  failed.  For  the  full  four-sensor  case,  the  cost  func- 

2 

tlon  has  minimum  value  of  6,38  EU  . With  one  sensor  failed,  this 

2 

value  increased  to  10.63  EU  . 


The  results  of  aligning  one  sensor  spin  axis  along  x and  the  three 

P 

others  in  an  umbrella  are  given  in  Figure  10,  For  the  four-operational 

2 

sensor  case,  the  cost  function  had  minimum*  value  equal  to  6.4  EU  . 

With  one  of  the  umbrella  sensors  failed,  the  minimum  value  of  J in- 
2 

creased  to  10.4  EU  at  an  umbrella  angle  of  15  deg.  It  is  interesting 
to  note  that  the  cost  function  vs  umbrella  angle  curve  of  case  (2)  re- 
mains relatively  flat  over  a fairly  large  range  of  umbrella  angles. 

This  situation  would  allow  for  larger  instrument  installation  tolerances 
without  sacrifice  of  gradient  determination  accuracy.  Accurate  measure- 
ment of  the  spin  axis  alignment  would  be  required. 
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HUEHC5/BCLL  ERflDIQHCTER  PERFORMnNO: 


FIG 


HUGHE5/EEU.  ERflDlONETER  PERrORHIiNCC 


THREE  SPIN  AXES  IN  UMBRELLA,  FOURTH  SPIN  AXIS  ALONG  x . UMBRELLA 
SENSOT  SPIN  AXES  SYMMETRICALLY  PLACED  ABOUT  x (120  deg’^intervals) 


MEASURING  THE  GRAVITY  GRADIENT  TENSCm  WITH  THE 
DRAPER  LAB  GRAVITY  GRADIOMETER 


An  orientation  analysis  similar  to  that  carried  out  for  the  rotating 
gradiometers  was  also  carried  out  for  the  Draper  Lab's  static  instrument 
system.  For  a gradiometer  triad  composed  of  three  spherical  sensors, 
the  defining  equation  has  the  form  of  (47) , Matrix  M of  (47)  is  similar 
in  nature  to  M'  of  (28)  in  that  it  is  composed  of  products  of  Euler 

rotations.  Note  that  the  scalar  multiplicative  coefficient  for 

2 2 

this  instrument  is  2mr  , whereas  it  was  4mr  for  the  Hughes'  sensor. 


= 2mr 


l_l  1 1-0  0 oJ[_ 


The  analytical  development  of  the  gradient  tensor  estimate  for  this  instru- 
ment system  is  similar  to  that  of  the  Hughes'  instrument  and  appears  in  Ref.  6. 
Optimal  sensor  orientation  is  again  defined  to  be  that  which  minimizes 
the  cost  function  specified  in  (38).  The  results  of  the  three-instrument 
analysis  are  shown  in  Fig.  11.  For  the  Draper  Lab's  sensor,  the  between- 
the— weights  axis  has  essentially  the  same  significance  as  the  spin  axis  has 
for  the  rotating  sensors.  The  J axis  angular  rotations,  P,  are 
identical  to  the  umbrella  angles  discussed  in  connection  with  the  Hughes' 
sensor  spin  axis  orientations. 


-21- 


The  performance  plots  of  Fig.  11  show  that  the  cost  function  J is 
minimized  when  the  between-the-weights  axes  are  symmetrically  located 

✓V  As  ^ 

about  X and  rotated  out  of  the  y - z plane  through  an  umbrella 

p p p 

angle  of  approximately  58  degrees.  The  value  of  J at  that  point  was 

2 

found  to  be  approximately  10.8  EU  . Note  that  for  the  symmetric  between- 
the-weights  axes  case,  a singularity  exists  at  an  umbrella  angle  of 
35.2644  deg.  It  is  easy  to  show  that  a singularity  always  exists  for 
this  three  sensor  system  when  the  between-the-weights  axes  form  an  orthog- 
onal triad  [6]. 

Four  sensor  gradient  tensor  estimation  was  also  considered  for  the 
Draper  Lab's  sensor.  The  following  cases  were  considered: 

(1)  Four  sensors  symmetrically  located  about  x^  having  identical 
umbrella  angles  j5; 

(2)  A single  sensor  having  its  between-the-weights  axis  vertical. 

Three  sensors  symmetrically  located  about  having  identical 

umbrella  angles  0, 

The  results  of  these  tests  for  the  four— sensor  and  one-sensor-f ailed 
cases  are  given  in  Figs.  12  and  13. 

Tile  cost  function  vs  umbrella  angle  curve  for  the  case  where  all  four 

As 

between-the-weights  (BT\V)  axes  are  symmetrically  spaced  about  x^  appears 

in  Fig.  12.  With  four  sensors  operating,  the  minimum  cost  function  occurred 

2 

at  an  umbrella  angle  of  50  deg  and  had  a value  of  7,1  EU  . W'ith  one 

sensor  failed  , the  minimum  value  of  the  cost  function  increased  to  11.9 
2 

EU  at  55  deg. 

The  effect  of  removing  one  of  the  BTW  axes  from  the  umbrella  and 

placing  it  along  x^  is  shown  in  Fig.  13,  For  this  case,  the  minimum 

value  of  J occurs  at  an  umbrella  angle  of  20  deg  and  has  value  equal 

to  8,25  EU^.  At  3 equal  to  55  deg,  the  value  of  J is  8,6  EU^, 

With  one  umbrella  sensor  failed,  the  minimum  value  of  the  cost  function 

2 

was  found  to  be  12.66  EU  at  an  umbrella  angle  of  35  degrees. 
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FIG.  11  GRADIENT  MEASUREMENT  USING  THREE  DRAPER  LAB  GRADiaMETERS 


<SI>RL  BRAD  IOMETER  PERFORMANCE 


G.  CONCLUSIONS 


The  results  presented  in  this  paper  show  that  three  spinning  gravity 
gradlometers  will  provide  a gradient  estimate  having  minimum  error  when 
the  sensor  spin  axes  form  an  orthogonal  triad,  A fourth  gradiometer  pro- 
vides redundancy  and  prevents  mission  loss  due  to  a single  instrument  fail 
ure.  It  has  been  shown  that  with  four  spinning  sensors,  the  best  informa- 
tion is  obtained  with  an  umbrella  configuration  in  which  one  sensor 
spin  axis  lies  along  x , and  the  other  spin  axes  are  symmetrically 
spaced  about  x^.  It  was  shown  that  substantial  errors  (on  the  order  of 
5 deg)  in  instrument  mounting  will  not  substantially  degrade  the  perform- 
ance of  the  sensor  system.  Furthermore,  it  is  concluded  that  for  a 
spinning  sensor  system,  an  umbrella  angle  of  20  deg  provides  a reasonable 
compromise  between  four  sensor  system  optimal  performance  and  performance 
degradation  due  to  the  loss  of  one  sensor. 

For  the  three  static  gradiometer  system,  it  was  found  that  BTW  axis 

✓V 

symmetry  about  x^,  and  an  umbrella  angle  of  approximately  58  deg  mini- 
mized the  cost  function  J,  This  result  was  in  general  agreement  with 
those  of  Ref,  4,  It  is  also  concluded  that  a four-static  gradiometer 
system  be  configured  with  one  BTW  axis  along  x^  and  the  others  symmetri- 
cally located  about  x^  with  umbrella  angles  of  20  deg.  The  20  deg 
umbrella  angle  selection  is  based  on  one-instrument-failure  considerations 
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Chapter  II 
STRESS  RELAXATION 


The  phenoaena  of  creep,  creep  recovery,  relaxation,  microstrain,  or 
whatever  term  is  applied,  aje  described  in  a variety  of  ways,  the  reason 
for  which  is  apparent  disagreement  of  the  actual  mechanism  involved. 
Johnson  [7]  considers  metal  as  an  aggregate  of  crystalline  domains 
embedded  ina  viscous  medium  (the  intercrystalline  boundary).  Assuming, 
for  the  moment,  three  stages  of  creep,  the  first  stage  is  postulated  to 
occur  by  viscous  motion  of  the  Intercrystalline  boundary.  The  second 
stage  occurs  by  slipping  of  the  crystalline  domains  within  the  boundary, 
and  the  third  stage  by  fragmentation  of  the  crystalline  domains. 

As  material  grain  size  Increases,  creep  rate  decreases  due  to  the 
reduction  in  the  number  of  Intercrystalline  boundaries.  The  extension  of 
this,  of  course,  is  that  single  crystals  would  exhibit  no  viscous  creep. 

Orowan  [8]  claims  single  crystals  do  creep  due  to  slip  of  atoms 
within  the  crystalline  lattice.  He  further  asserts  that  the  stress  at 
which  a single  crystal  begins  to  yield  depends  on  the  amount  of  prior 
strain,  i.e.,  the  strain  history. 

Creep  is  commonly  described  as  occuring  in  three  stages: 

(1)  The  primary  stage  during  which  strain  decreases  rapidly  with 
time ; 

(2)  The  secondary  stage  when  the  strain  rate  remains  fairly  con- 
stant ; 

(3)  The  tertiary  stage  in  which  strain  Increases  rapidly. 

During  the  primary  stage,  the  creep  rate  decreases  rapidly  due  to 
a buildup  in  crystal  dislocations  which  pile  up  and  inhibit  the  viscous 
flow. 

The  secondary  stage  is  a tug-of-war  between  strain  hardening  and 
reduction  in  specimen  cross  section  resulting  in  an  essentially  constant 
strain  rate  [ 9 ] . 
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For  analysis  purposes,  materials  that  creep  are  represented  as  one 
of  three  models: 


Maxwell  model  Is  represented  by  a linear  spring  and  linear  dash- 
pot  connected  In  series,  vnien  loaded,  such  a material  responds  with 
Instantaneous  elastic  strain  followed  by  continuously  creeping 
strain.  Unloading  results  In  only  partial  recovery  of  strain. 


Kelvin  model  Is  represented  by  a linear  spring  and  linear  dash- 
pot  connected  In  parallel.  A Kelvin  material,  when  loaded,  first 
begins  to  creep  at  a rate  that  decreases  to  zero.  Unloading 
results  In  gradual  recovery  to  the  original  size. 


Three  element  or  Burgers  model  Is  simply  a Maxwell  and  a Kelvin 
model  connected-on  series.  This  type  material  displays  instantaneous 
elastic  response,  creep  with  time,  elastic  memory;  and  permanent 
deformation  under  load.  Removal  of  the  load  results  In  partial 
strain  recovery. 

Creep  recovery  Is  most  easily  explained  In  terms  of  viscous  be- 
havior. For  most  polycrystalllne  metals,  the  amount  of  creep  recovery 
Is  very  small,  whereas  In  plastics.  It  Is  more  pronounced. 

The  time  required  for  the  load,  when  held  at  constant  strain,  to 
relax  to  l/e  of  Its  original  value  Is  termed  the  'relaxation*  time 
[18],  [19].  Thus  a viscoelastic  material  Is  characterized  by  moduli 
that  are  time  dependent  and  type-of**loadlng  dependent.  Since  most  visco- 
elastic materials  behave  linearly  over  some  range.  It  Is  possible  to 
Identify  those  limits  In  terms  of  stress,  strain,  and  time,  and  the 
8 tress-strain- time  behavior  of  the  material  can  be  fairly  accurately 
predicted  for  design  purposes  [20]. 

Materials  for  use  In  gravity  gradlometer  design,  regardless  of 
what  they  are,  will  exhibit  some  degree  of  creep,  creep  recovery,  or 
stress  relaxation  and  this  effect  must  be  taken  Into  account  [21]. 
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